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Abstract. High energy electrons and positrons from decaying dark matter can 
produce a significant flux of gamma rays by inverse Compton off low energy photons 
in the interstellar radiation field. This possibility is inevitably related with the dark 
matter interpretation of the observed PAMELA and FERMI excesses. The aim of this 
paper is providing a simple and universal method to constraint dark matter models 
which produce electrons and positrons in their decay by using the Fermi LAT gamma- 
ray observations in the energy range between 0.5 GeV and 300 GeV. We provide a set 
of universal response functions that, once convolved with a specific dark matter model 
produce the desired constraint. Our response functions contain all the astrophysical 
inputs such as the electron propagation in the galaxy, the dark matter profile, the 
gamma-ray fluxes of known origin, and the Fermi LAT data. We study the uncertainties 
in the determination of the response functions and apply them to place constraints on 
some specific dark matter decay models that can well fit the positron and electron 
fluxes observed by PAMELA and Fermi LAT. To this end we also take into account 
prompt radiation from the dark matter decay. We find that with the available data 
decaying dark matter cannot be excluded as source of the PAMELA positron excess. 
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1. Introduction 

The existence of dark matter in our Universe is nowadays a widely accepted fact. Dark 
matter constitutes the dominant fraction of matter present in the Universe. However, 
all the astrophysical evidences for dark matter, such as those coming from gravitational 
lensing, galaxy rotation curves and cosmic microwave background (CMB) anisotropies, 
are purely gravitational [1], and its particle nature remains unknown. 

The most popular type of dark matter (DM) candidate, a weakly interacting 
massive particle (WIMP), can naturally reproduce the observed DM abundance due 
to effective self-annihilation in the early Universe. Today, this same annihilation 
process could produce an observable contribution to the measured cosmic-ray fluxes 
on Earth [2, 3, 4, 5] and have impact on cosmological observations [6]. However, 
such an indirect detection of DM is also possible if DM decays at a sufficiently large 
rate. Indeed, there exists a number of interesting and theoretically well motivated 
DM models that predict the decay of DM on cosmological time scales, namely with 
lifetimes around and above r x ~ (9(10 26 s), which are typically required to be not 
in conflict with current observational limits [7, 8, 9]. Among these candidates is the 
gravitino in small -R-parity breaking scenarios (motivated by requiring a consistent 
thermal history of the Universe [10]) and models of sterile neutrinos, whose long lifetime 
is due to tiny Yukawa couplings (see Ref. [11] and references therein). Other interesting 
models include kinetically mixed hidden gauge bosons and gauginos [12, 13]. Even 
in models where the DM candidate is stable in the first place, the consideration of 
higher-dimensional operators often renders the DM candidate particle unstable with 
cosmological lifetimes [14]. Since the indirect detection signals from decay differ in 
general from the ones of annihilation, a dedicated study of decaying DM signals is 
warranted . 

On the observational side, recent results by several Cosmic Ray (CR) and 7- 
ray experiments have seriously challenged the standard description and interpretation 
of the generation of high-energy CR electrons, positrons and 7-rays, thus giving 
room also to DM interpretations. In particular, the PAMELA satellite reported the 
presence of an increase in the positron fraction J(e + )/J(e + + e~) above 10 GeV [15]. 
If there is no primary source of CR positrons in the Galaxy, and since secondary 
positrons originate from the interactions between CR nuclei (mainly protons and 
Helium) and the interstellar gas, their spectrum on Earth should be softer than the 
primary electron spectrum. Therefore, the positron fraction is expected to decrease 
with increasing energy. The observed rise at high energy (the so called positron 
fraction excess) leads then to postulate the existence of an additional primary source of 
positrons with energies ranging from 10 to at least 100 GeV. Although considering 
more carefully standard astrophysical sources, such as pulsars, provides a possible 
explanation for this excess [16, 17, 18, 19, 20], also dark matter interpretations are 
possible. In particular, DM decays might provide the required amount of primary 
positrons [14, 12, 21, 22, 23, 24]. However, we should stress that the requirement of 
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an additional primary positron source is not strictly necessary to interpret the data. 
The positron excess can be interpreted also as due to propagation effects, either in the 
standard sources or in the ISM. For example, the excess could be due to inhomogeneities 
of CR sources [25], or it could be a consequence of reacceleration of secondary cosmic rays 
in supernova remnants [26]. Also an increase of starlight and interstellar gas densities 
could reproduce the positron excess [27]. Lastly, also nested leaky box models can in 
principle provide a good fit to the positron fraction [28, 29]. 

Apart from the positron fraction, also the e~ + e + spectra from the Fermi LAT 
measurement [30] and HESS [31] are intriguing. The Fermi LAT instrument measured 
a total e + + e~ spectrum at high energies from 20 GeV to 1 TeV on Earth significantly 
harder than expected according to propagation models, and with a hint of the presence 
of a broad bump at a few hundred GeV, which is smoothly connected with HESS 
observations around the TeV (see, e.g. [32, 33] for a review of possible interpretations of 
these observations). This also provides a motivation to explain the excess using decaying 
dark matter which can easily produce large numbers of energetic positrons and electrons. 
Both these striking observations have motivated a plethora of works on both building 
and constraining new dark matter models . 

High energy CR electrons and positrons produce 7-rays through the processes of 
inverse Compton scattering (ICS) off low energy photons in the interstellar radiation field 
(ISRF) and through bremsstrahlung emission due to the interaction with the interstellar 
medium (ISM). The Fermi Gamma- Ray Space Telescope is currently observing the 7- 
ray sky in the energy range between 30 MeV and 300 GeV with unprecedented angular 
resolution and precise energy sensitivity. Data obtained during the first year of operation 
of Fermi are about to be released officially, and preliminary data on the 7-ray spectrum 
in some regions of the sky were already made public in several Conferences. It is then 
timely to provide tools to exploit this detailed information in an efficient way. 

At present, most attempts to construct dark matter candidates as the interpretation 
of the PAMELA and 7-ray data are inevitably model dependent. Thus, it is useful 
to develop model independent analysis tools such as response functions. These are 
functions of injected electron or positron energy which are obtained from comparing 
the 7-ray spectrum resulting on Earth after propagation (which we call signal) with the 
observed 7-ray spectrum (which henceforth will be called background) . Indeed, this 
method has already been used to constrain decaying DM scenarios with the observed 
galactic synchrotron radiation maps [34]. Since the propagation equation is linear in 
the source term, each injected electron energy evolves independently. Therefore, with a 
finite number of numerical simulations at different injection energies we can construct 
a numerical e ± -response function of signal-to-background . The e ± -response function 
encodes our knowledge about the propagation of CR electrons in the Galaxy, as well as 
about the distribution of DM, but is independent of the particle physics model leading 
to the actual spectra of electrons and positrons after DM decay. Moreover, the e ± - 
response functions we compute in this work are normalized to the present 7-ray data. 
Hence, constraints can be obtained by the convolution of the e ± -response functions with 
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a given DM decay spectrum and requiring that the result be smaller than the product 
of DM mass and lifetime in suitable units. If this condition is violated, the DM model 
predicts 7-ray fluxes (from ICS and bremsstrahlung) that are higher than the data by 
at least 2a. 

In addition to 7-rays produced by ICS and bremsstrahlung of pairs resulting from 
DM decay, the prompt radiation of 7-rays during the DM decay can be a very relevant 
signature of DM decays . However, since prompt radiation depends on the particle 
physics, we cannot include it in our e ± -response functions and has to be included by 
hand when pursuing constraints on a specific model. 

The paper is organized as follows. In Sec. 2, we describe our method of e ± -response 
functions to provide constrains on any dark matter model involving electrons and 
positrons in its decay products. In Sect. 3 we present the e ± -response functions based 
on the Fermi data and discuss how they are improved by the removal of astrophysical 
contributions of known origin to the 7— ray signal. In Sect. 4 we apply our e ± -response 
functions to some specific dark matter decay scenarios and compare it to the contribution 
of prompt radiation to the 7— ray constraints. Finally, in Set. 5 we present our final 
discussions and summarize. 



2. Gamma-rays from Decaying Dark Matter 

In this section we will discuss the different sources of gamma radiation resulting from 
dark matter decay. Generically, 7-rays are expected as a secondary product of the 
interactions of CRs with the galactic gas and radiation field. It is well known that 
CRs in the energy range we are interested in obey the following propagation equation 
(Ginzburg & Syrovatskii [35]) 

9N e± „ XAre ±, d (. p„ \ Are± d _ 2 d N e± 



at -v-a>v-«.)jv -^*ir = 

= Q e± (p,r, z) + c£n g „(r, z) c r ± N CR . (1) 

Here N e± (p, r, z) is the number density of the electrons/positrons; p is their momentum; 
j3 their velocity in units of the speed of light c; <r e ± is the production cross section of 
electrons/positrons from spallation of CR protons and heavier nuclei with density N CR 
onto the ISM gas, whose density is n gas ; D = (3D (1Z/GV) S is the spatial diffusion 
coefficient which is a function of the particle rigidity TZ; v c is the convection velocity; 
p are continuous energy losses, whose relevant timescales r = p/p are shown in Fig. 1. 
The last term on the l.h.s. of Eq. (1) describes diffusive reacceleration of CR in the 
turbulent galactic magnetic field. In the quasi-liner theory the diffusion coefficient in 
momentum space D pp is related to the spatial diffusion coefficient by the relationship 
(see e.g. [36]) 

D - A - ?A]? (2) 

pp ~ 35(4 -52) ( 4 -6) D [ ) 
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Figure 1. Energy loss time scales for CR electrons at Solar position (left panel) and 
at r = and z = 2kpc (right panel). The black solid curve represents IC losses, cyan 
long dashed is for synchrotron losses, red short dashed for bremsstrahlung losses, blue 
dotted for Coulomb losses, green dot-dashed for ionization losses. 

where va is the Alfven velocity. Finally, Q e± (p,r, z) represents the primary injection 
term. 

The electron distribution can be computed numerically in the stationary limit 
dN e /dt = by using the methods of [41, 42, 43]. However, for this paper we have 
used our own propagation codes, which have strongly improved convergence behavior 
for mono-energetic injection spectra. From this distribution we will then compute the 
7-ray spectra and maps produced by ICS and bremsstrahlung. 

We will firstly discuss inverse Compton scattering and bremsstrahlung radiation, 
and secondly prompt radiation, both with their galactic and extragalactic contributions, 
with particular emphasis on the related uncertainties. 

2.1. Inverse Compton scattering and bremsstrahlung 

Electrons and positrons produced in the decay of dark matter give rise to a 7-ray 
signal, coming from the inverse Compton scattering (ICS) of low energy photons of the 
interstellar radiation field (ISRF). A further, mostly subdominant, contribution to the 
7-ray flux comes from bremsstrahlung of the electrons and positrons when scattering 
with the galactic gas. In the context of the PAMELA positron excess, ICS radiation from 
decaying dark matter has been discussed in Refs. [37, 38], for the case of annihilating 
dark matter see Refs. [39, 40]. 

The ICS 7-ray flux due to injected mono-energetic electrons and positrons from 
dark matter decay inside of our Galaxy, and arriving on Earth from the direction Cl 
(characterized by galactic latitude b and longitude I), is given by 

J lcs (Cl,E 7 ;E e ) = ^- I ols I dEP(E 11 E)n e (r,E;E e ) , (3) 

were n e (r, E; E e ) satisfies the stationary solutions to the propagation Eq. (1) for a 
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source term with monochromatic injection of positrons and electrons at an energy E e , 
P(Ej, E) is the differential photon emissivity and the spatial integral is along the line 
of sight. The relevant source function is given by 

Q(r,E) = ^-8(E-E e ). (4) 
r x m x 

where m x is the DM particle mass and r x is its lifetime. 

As dark matter profile p x (r) we take (assuming, as customary, spherical symmetry) 
the NFW profile [44] 

which we normalize to p© = 0.36 GeV/ cm 3 at the position of the Sun, r = 8.5 kpc, 
and we adopt r s = 20 kpc. Uncertainties from the halo profile will be discussed below. 
The differential emissivity P(E 1 , E) corresponds to two processes, 

P(£ 7 , E) = P 1C (E 7 , E) + P brcmss (£ 7 , E) . (6) 

The first term Pjc{E 1 , E) corresponds to inverse Compton scattering, which is derived by 
convolving the differential number density of target photons with the well known Klein- 
Nishina cross section. The second term, P bremss (i? 7 , E), stems from the bremsstrahlung 
emission due to deflection of relativistic electrons and positrons in the electrostatic 
potential of interstellar gas atoms and molecules. Since this contribution is subdominant 
in the energy range we consider in this work, we will always refer to the ICS channel as 
our reference channel in the following. More details of the calculation can be found in 
[45]. 

The galactic radiation field in the galactic plane consists mainly of starlight at 
wavelength A ~ 1 /im, diffuse dust emission at A ~ 100 /im and the CMB. Notice 
however that since most of the radiation (CMB aside) is produced in the galactic plane, 
the actual radiation fields at large latitudes are poorly known, and the computation of 
ICS radiation from high latitudes suffers from large uncertainties. In this work we will 
adopt the ISRF presented in Ref. [47]. 

The predicted 7-ray flux crucially depends on the propagation model. We list 
in Tab. 1 three different models characterized by different choices of the propagation 
parameters in Eq. (1). By comparing the results obtained with these three models 
we will estimate the uncertainties associated to the poor knowledge of the propagation 
models. The MIN and MAX models are diffusion-reacceleration-convection models. They 
were proposed in [48, 49, 50] and are known to be compatible with the observed B/C ratio 
and produce the minimum and maximum observationally allowed fluxes of antiprotons 
on Earth, respectively. The LI model is derived from [51], in which an energy dependent 
analysis of recent data about secondary /primary ratios allowed a fairly accurate study 
of the diffusion parameters 5 and Dq/L.\ Although the main aim of [51] was not to find 
a best fit model able to reproduce all the observed spectra, it is remarkable that the 

% It is a general fact that stable secondary/primary ratios do not allow to probe separately the 
magnitude of the diffusion coefficient and the height of the diffusion region. Unstable/stable ratios, 
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km/s/kpc 
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MIN 


0.85/0.85 


0.048 


20 


1 


13.5 





22.4 


0.1 


LI 


0.50/0.50 


4.6 


20 


4 








10 


4 


MAX 


0.46/0.46 


2.31 


20 


15 


5 





117.6 


0.1 



Table 1. Typical combinations of diffusion parameters that are consistent with an 
analysis of CR nuclei secondary /primary ratios. The MIN and MAX propagation 
models correspond to minimal and maximal primary antiproton fluxes, respectively, 
while the LI model can provide a good description of B/C, p/p and data on other 
secondary/primary ratios above 1 GeV/n. 

1 Below/above the break in rigidity at K = 4 GV for the MIN, MAX and LI model. 

diffusion parameters determined via an high energy analysis are able to describe data 
down to energies of the order of 1 GeV/nucleon , and also to reproduce with reasonable 
accuracy the PAMELA measurements of the antiproton flux . It must be noticed, 
however, that this is only one out of many different possible models in agreement with 
nuclei CR observations (also simple leaky-box models can in fact succeed in reproducing 
high energy data [52]). In fact, all three models listed in Tab. 1 well reproduce available 
data, in spite of the different parameters they adopt. Moreover, all these models assume 
cylindrical symmetry and diffusion to be homogeneous and uniform over the whole 
Galaxy, which is a rough approximation of the galactic propagation regime. Yet, this 
simple modeling allows to reasonably reproduce available data with relatively few free 
parameters. A study of the effects of a possible spatial dependence of the diffusion 
coefficient in a general three-dimensional geometry is beyond the scope of this work. 

For illustrative purposes, we show in Fig. 2 the predicted 7-ray emission for the 
LI model of Tab. 1 at 0.1, 1 and 10 GeV, respectively, assuming a dark matter decay 
into e + e~, where m x = 200 GeV and r x = 10 26 s. We also show for comparison the ICS 
radiation from primary electrons of astrophysical origin. In general, as apparent from 
these plots, dark matter induced ICS radiation extends to higher latitudes than the 
ICS radiation due to the contribution of standard astrophysical source, which is mainly 
concentrated on the galactic disk. 

In addition to the ICS radiation produced in our Galaxy there is also a related 
extragalactic contribution, resulting from scattering of electrons from dark matter 
decaying outside of our Galaxy with the CMB. This component can potentially dominate 
the overall ICS fluxes at low energies, and we include it for completeness. The calculation 
is straightforward and details can be found in, e.g., Ref. [37, 40]. In contrast to 
these references, we also took into account absorption effects due to inelastic scattering 
between ICS photons and the intergalactic background light (IBL), following Ref. [53, 38] 
(adopting the "fast evolution" model). 

such as the 10 Bc/ 9 Be, can in principle provide such a discrimination. The available data on such ratios, 
however, have very large errors, so that it is extremely difficult to extract information from them. 
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Figure 2. The 7-ray emissions at 0.1 GeV, 1 GeV and 10 GeV (upper panel from 
left to right) produced by dark matter particles decaying into e + e~ pairs, where 
m x = 200 GeV, t x = 10 26 s. Results hold for the LI diffusion model of Tab. 1 and for 
the NFW halo profile. The lower panel shows the ICS radiation from astrophysical 
sources at 10 GeV for comparison (again from model LI). The color scaling corresponds 
to the logarithm to the base 10 of the flux in GeV/s/cm 2 /sr. Note that the color scale 
corresponds to the same flux range in all panels. 

Another contribution to the ICS radiation comes from electrons and positrons 
produced in the dark matter decay inside of our dark matter halo, but outside of the 
diffusion zone. If these particles are far enough from the galactic disk region, beyond a 
few kpc, their main energy loss channel is scattering with the CMB. The corresponding 
propagation length is (9(100 kpc) for 1 TeV particles, and O(10 kpc) for 10 TeV particles. 
Calculating the corresponding ICS flux while neglecting the motion of the particles in 
general overestimates this flux by an 0(1) factor. Being conservative we do not include 
this radiation component to our bounds and leave a more careful calculation for future 
work. 

Uncertainties. The largest uncertainties in the predicted 7-ray fluxes come from poorly 
known propagation parameters, in particular from the height of the diffusion zone. The 
corresponding uncertainties can reach one order of magnitude. 

As we already noticed, the height L is only poorly known since its determination 
is degenerate with the diffusion coefficient. The most widely adopted range of variation 
of L is between 1 and 15 kpc, based on the poor quality data on 10 Be/ 9 Be. Due to this 
uncertainty it is hard to obtain any definitive constraint on, e.g., the supersymmetric 
parameter space of DM based on current anti-proton data [49, 54]. The forthcoming 
AMS-02 experiment [55] will however provide very accurate data on unstable/stable 
ratios (as well as for B/C and other stable secondary /primary ratios), which might 
allow a more precise determination of the diffusion height scale L. 
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In our case it is clear that a larger height of the diffusion zone leads to more DM 
decays contributing to the 7-ray flux, because more electrons of DM origin are confined 
in the diffusive region. Uncertainties on other CR propagation parameters, such as 
the Alfven velocity v& and the convection velocity v c , are less relevant in affecting the 
electron distribution of DM, and especially affect only the electrons below around 10 
GeV since higher energy electrons lose energy too rapidly via ICS and can not propagate 
over long distances. 

We found that the full-sky ICS 7— ray emission induced by very high energy 
electrons from DM decay obtained in the LI model is comparable with the one obtained 
using other widely known models, namely DC and DR [41] which adopt the same height 
of the diffusion zone as our LI model. 

Another source of uncertainty comes from the dependence of the 7-ray emission 
on the halo profiles. We compared the fluxes predicted for different halo models (with 
parameters as in Ref. [38]) and found that for shallower halo profiles such as the Kra [56] 
and the isothermal [57] profile the 7-ray emission is reduced by around 10%. On the 
other hand, in case of the Einasto profile (see [3] and references therein) the flux 
is enhanced by 30%. Since this variation is subdominant when compared with the 
uncertainties of the propagation models, we will simply adopt the NFW halo profile in 
the rest of the present work. 

The 7-ray emission also depends on the Galactic magnetic field, since synchrotron 
losses can be of the order of ICS radiation losses at high electron energies. The magnetic 
field profile close to the Galactic center is quite uncertain and could be considerably 
higher than a few /xG [59]. In this paper, we adopt the field from Ref. [58] which 
matches the 408 MHz synchrotron distribution. We compared our results with another 
widely used magnetic field and found that the 7-ray fluxes at high energies increase by 
just 15% by changing the magnetic field model to the one presented in Ref. [60] which 
is based on a large-scale data set on starlight polarization. 

2.2. Prompt Radiation 

Although this work focuses on ICS radiation produced by electrons and positrons from 
decaying dark matter, we cannot overlook the fact that the prompt radiation of gamma 
rays can be a very competitive signature of dark matter decays. Indeed, in many realistic 
cases it turns out that this component can be larger than the ICS and therefore it can 
actually give the strongest constraints. 

Prompt radiation from dark matter decay is produced more frequently inside our 
overdense Galactic dark matter halo§ but it is also produced at cosmological distances. 
At energies around 10 GeV or below, the magnitude of the halo and extragalactic fluxes 
are of the same order when looking in direction of the anti-galactic center, whereas at 

§ Although the halo profile is expected to be approximately isotropic, the corresponding flux at Earth 
exhibits a strong dipolc-likc anisotropy due to the offset between the Sun and the galactic center, which 
can be used to distinguish it from the extragalactic 7-ray background, see Ref. [38]. 
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higher energies around and above 100 GeV the inelastic scattering between 7 rays and 
the IBL reduces the extragalactic component significantly and can not be neglected. 

We include the galactic and the extragalactic prompt radiation, following the 
calculations outlined in Ref. [38]. Uncertainties come mainly from the adopted dark 
matter halo profile and its normalization at position of the Sun. 



3. e -response Functions from Inverse Compton Scattering 

In this section we derive the e ± -response functions for the 7-ray emission induced by 
DM decays into electrons and positrons by comparing our predictions with the 1-year 
observations of Fermi LAT in the energy range of 0.5 GeV to 300 GeV. 

As data we use the Fermi LAT 7-ray maps as derived in Ref. [61]. || In this analysis 
the events were binned into the eight energy ranges 0.5-1 GeV, 1-2 GeV, 2-5 GeV, 
5 - 10 GeV, 10 - 20 GeV, 20 - 50 GeV, 50 - 100 GeV and 100 - 300 GeV. 

A few words of caution are required concerning the data basis: The adopted 7-ray 
maps are based on the "diffuse" event class, which at energies above 50 or 100 GeV 
suffers from background contamination [62] that becomes relevant when diffuse fluxes 
are studied. In the highest energy regime the background contamination might be on 
the 50 - 80% level. Furthermore, no attempt of subtracting point sources were made. 
The adopted statistical errors are always derived from the exposure 3 x 10 10 cm 2 s, which 
is good enough for the purpose of this paper. ^[ We do not attempt to include systematic 
errors in the analysis, which in light of the large background contamination of the data 
would be too premature. 



3.1. Definitions and Optimal Sky Patch 

Following the same Ansatz we made for the case of synchrotron radiation proposed in 
Ref. [34], we introduce response functions, which are functions of the electron injection 
energy and are associated to 7-ray observations in a sky patch Afl and in an energy 
band E Q < E 1 < E\. The e ± -response functions are defined as the ratio of the predicted 
7-ray fluxes resulting from decaying dark matter to the observed fluxes as 

F £o*r A n. F \ = S. ^ dQJlCS ^ glj Ee ^ ( T x \ ( m X W 
7 ^^>~ J o ^ (AO)+2 . 5Jobs vl(F s Al00GeW Uj 

where J ICS (f2, E 7 ; E e ) is calculated for the injection spectrum Eq. (4) and we adopt the 
conservative attitude of adding, in each energy and angular bin, to the central value of 
the observed flux the associated 2a error. These functions depend on neither t x nor m x 

I We also performed an analysis of the publicly available event data on http://fermi.gsfc.nasa. 
gov/ssc/data/, coming to identical results as Ref. [61]. We decided to use the maps from [61] in order 
to keep the data basis of our analysis easily accessible. 

% We cross-checked with our own analysis of the publicly available event data that this gives indeed 
the correct number counts at high energies. 
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of Eq. (4), and constraints on a given DM decay model can then be easily cast in the 
form 

£^^ 1(Afi;Be) ^<(^)(_^_), (8) 

where dN e /dE e is the electron/positron spectrum obtained from DM decay in some 
specific particle physics model. We stress again that the advantage of the e ± -response 
function approach is that the e ± -response functions are independent of the specific 
DM decay spectrum. The method is hence directly applicable to any DM model, and, 
moreover, allows a discussion of the typical characteristics of ICS radiation from DM 
decay in a model-independent way. 

The e ± -response functions depend crucially on the chosen patch Af2, which in the 
ideal case should cover the area with the largest signal-to-background ratio, maximizing 
the e ± -response functions. For the sake of clarity of our approach, and to avoid problems 
with statistical bias related to adaptive methods (which are related to downward 
statistical fluctuations and become severe if the statistics is low, see Ref. [63]), we 
choose a fixed patch with large signal-to-background ratio by inspection of the signal- 
to-background maps. 

Clearly, the optimal region might in principle depend on the observed energy range 
and injection energy of electrons and positrons. But it turns out that for injection 
energies around 100 GeV - 10 TeV and high enough 7-ray energies, the optimal region 
is always located south of the galactic center. This situation is illustrated in Fig. 3, 
which shows an example of a signal-to-background map for observed energies between 
0.5 and 1 GeV, where the statistics is very good, for dark matter decay producing 
monochromatic electrons and positrons at 100 GeV. The map already suggests that the 
location of the optimal patch for constraining inverse Compton light from decaying dark 
matter actually lies in a region close to the galactic center, located south of the galactic 
plane. 

Before discussing this further we note that an exception occurs for very high 
injection energies in the 1-10 TeV region, and low enough observed 7-ray energies. 
There the overall ICS flux can actually be dominated by extragalactic ICS contributions 
from scattering between electrons and positrons from dark matter decay with the CMB. 
In these cases the optimal patch would be located at the pole regions. The same holds in 
general true for prompt radiation from dark matter decay in and outside of the Galaxy, 
which has a much shallower angular profile than the galactic ICS component. 

A more quantitative description of the situation can be found in the plots shown 
in Fig. 4. The black lines in the upper four panels show the signal-to-background ratio 
in different observed energy ranges as function of galactic latitude or longitude. The 
injection energy is now fixed to 1 TeV, but the results stay qualitatively the same for 
other injection energies. As expected from the above discussion, at high energies {upper 
panels) the signal-to-background ratio is maximal in a region close to the galactic center, 
whereas at very low energies {middle panels) it is maximal at high latitudes, due to the 
extragalactic ICS component. The same dominance at high latitudes is also present in 
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Figure 3. Signal-to-background map of ICS radiation from dark matter with 
to x = 200 GeV, r x = 10 26 s decaying into e ± -pairs, compared to the Fermi LAT 
7-ray observations in the 0.5 - 1 GeV regime. Results hold for the LI propagation 
model of Tab. 1. Note the logarithmic color scaling, warmer colors indicate larger 
signal-to-background . 



case of prompt radiation (lower panels). 

For the derivation of the e ± -response functions we will concentrate on the patch S 
close to the galactic center and defined by |/| < 20° and —18° < b < —10°, which is 
marked by the colored region in Fig. 4. We checked that this patch indeed maximizes 
the obtained constraints when varying the patch boundaries, except for the very highest 
energy region (100 GeV - 300 GeV), where however the statistical error is large. For 
this patch, we construct the e ± -response functions by performing simulations injecting 
different mono-energetic electrons and comparing the results with observations according 
to Eq. (7). 

Given that in some cases the optimal patch is actually located at the galactic 
pole regions we will also calculate and present bounds that come from comparing 
the preliminary extragalactic 7-ray background as determined in Ref. [62] with the 
extragalactic ICS and isotropic prompt radiation component, see below. 

3.2. Response Functions without Foreground Subtraction 

Our results for the e ± -response function are shown in Fig. 5 as function of the 
electron/positron injection energy, for the 8 different energy ranges of Fermi LAT 
skymaps from Ref. [61]. The highest energy range provides the strongest constraint 
on decaying dark matter with very hard electron/positron energy spectrum. However, 
for lower injection energies in the 100 GeV— 1 TeV region, several energy ranges give 
actually roughly the same constraints. 

To illustrate the large uncertainties related to inverse Compton radiation from dark 
matter decay inside the diffusive halo, we show in Fig. 6 the e ± -response functions based 
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Figure 4. Signal-to-background ratios as function of galactic latitude (left panels) and 
longitude (right panels). The upper and middle panels correspond to pure ICS signal, 
the lower panels correspond to the pure prompt signal for comparison. Extragalactic 
and galactic radiation are taken into account. The black lines take into account 
as the whole observed signal, the green lines are obtained after subtraction of our 
reference model for the astrophysical component (Model LI). We find that the signal- 
to-background ratio of ICS radiation at higher 7-ray energies is maximized in the 
region <S defined by |Z| < 20° and -18° < b < -10°, which is indicated by the light 
red shaded region. 



on the highest and lowest 7-ray energy ranges for our three reference propagation models 
from Tab. 1. As emphasized before, the uncertainties on the e ± -response functions 
are dominated by the propagation model, especially for the lower energies, below 10 
GeV injection energy, where also effects of reacceleration become relevant. For higher 
injection energies above 10-100 GeV, where the response functions become of 0(1) and 
are hence relevant for the actual bounds, the uncertainties mainly stem from the height 
of the diffusion zone. In effect, high energy electrons and positrons lose energy in a 
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Figure 5. The e ± -response function F 7 based on 7-ray emission for the LI model 
of Tab. 1. The e ± -response functions are derived from the eight 7— ray energy ranges 
0.5 - 1 GcV, 1 - 2 GcV, 2-5 GeV, 5-10 GcV, 10 - 20 GeV, 20 - 50 GcV, 50 - 100 
GeV, and 100 — 300 GeV from top to bottom at left side, respectively. The underlying 
sky patch S is defined by |Z| < 20° and -18° < b < -10°. 




E(GeV) 

Figure 6. The propagation model dependence of the e ± -response function F 1 based on 
our fixed patch for the 7— ray energy range 0.5 — 1 GeV {yellow band, curves extending 
to low energies) and 100 — 300 GeV (black band, curves cutting off around 100 GeV). 
The width of the bands represents the variation within the MIN (green), LI (red) and 
MAX (blue) propagation models of Tab. 1. 
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Figure 7. The dependence of the e ± -response functions on subtraction of astrophysical 
contributions to the 7— ray signal. The e -response function for the LI model based 
on the raw observed map (solid) and on residual maps with 7— rays from ir° decay 
(dotted) and from all astrophysical processes (dashed, see text) removed. Red lines 
extending below 1 GeV are based on 7— ray flux observed in the energy range 0.5 — 1 
GeV and green lines are based on the interval 100 — 300 GeV. 

very short time compared to the diffusion time, thus making the other details of the 
propagation irrelevant. The MAX propagation model gives the strongest constraints due 
to its large diffusive halo, whereas the MIN propagation model minimizes the constraints. 
Moreover, for the MAX model re-acceleration shifts lower energy electrons to higher 
energies. This effect is however only relevant for electrons below around 10 GeV, and 
thus increases the 7-ray emission only in the MeV regime. Note that for the highest 
observed 7-ray energy region (100-300 GeV), one clearly finds a sharp cut off at low 
injection energies since 7— rays at such high energies cannot be produced from ICS of 
electrons/positrons injected at energies lower than 100 GeV. 

3.3. Response Functions with Subtraction of Astrophysical Foregrounds 

The e ± -response functions discussed so far are conservative because we did not 
attempt to subtract any astrophysical contribution to the 7-ray flux. In order to 
understand the conventional astrophysical 7-ray flux one needs to estimate the 7-ray 
emission from different galactic components. The most relevant production channels 
are nucleus-nucleus (mainly proton-proton) photoproduction via tt° decay and ICS 
and bremsstrahlung of CR electrons and positrons. It is generally found [58] that 
hadronically generated 7-rays dominate the flux at energies between 0.1 and 100 GeV 
and in the vicinity of the galactic plane, where most of the interstellar gas is located, 
while at lower and higher energies and at high latitudes ICS becomes comparable and 
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can dominate. Bremsstrahlung is usually a subdominant component. 

In Fig. 7 we show the e ± -response function for the LI model based on residual 7-ray 
maps obtained by subtracting 7-rays produced via tt° decay ICS and bremsstrahlung. 
In this foreground model the electron flux is adjusted to always lie below the electron 
flux observed by Fermi LAT, with a spectral index of around -3.2. The subtraction 
affects the results at low energies. For example, at E e ~ 10 TeV the e ± -response 
functions based on 7— ray fluxes observed at energies 0.5 — 1 GeV are increased by 
a factor of around five when subtracting the total astrophysical "foreground" at these 
7— ray energies. Again at E e ~ 10 TeV the e ± -response functions based on 7— ray fluxes 
observed at energies 100 — 300 GeV are increased by < 10% and ~ 15% by the removal 
of 7-rays originating from tt° decay and from all astrophysical processes, respectively. 
This demonstrates explicitly that constraints on dark matter decay can be improved 
by taking into account the removal of astrophysical contributions mentioned above. At 
high energies the removal turns out to be quite insufficient which is at least in part 
related to the strong background contamination in the adopted data. This situation 
will be improved when data with better background rejection becomes available. 

4. Constraint on Dark Matter Models 

In the previous sections we constructed e ± -response functions for DM decays which are 
independent of the particles physics details of the decay and only depend on the spatial 
distribution of the DM particles and the propagation of the produced electrons and 
positrons. These e ± -response functions are very useful to estimate the constraints on 
DM decay based on 7-rays produced by ICS within some specific DM model. However, 
within a given DM model 7-rays can be produced not only as secondaries of electron and 
positron propagation in the Galaxy, but also as final-state, or prompt, radiation arising 
in the decay [64, 65]. For selected decay channels, we calculated the corresponding 
galactic and extragalactic prompt fluxes in our selected patch as described above and 
added them to the one described by the e ± -response functions in order to derive 
constraints that come from the total prompt + ICS radiation flux of different dark 
matter decaying models. Other works on constraints on the decaying/ annihilating dark 
matter interpretation of the PAMELA positrons excess with recent Fermi LAT 7-data 
can be found in Ref. [24, 63, 66]. 

To illustrate the interplay of prompt and ICS radiation bounds, we show them in 
detail for the exemplary decay channel into in Fig. 8. In this plot, each of the 

green lines corresponds to bounds coming purely from the galactic and extragalactic ICS 
radiation (calculated from our response functions) for the eight different energy regimes 
of the data. Regions below the green lines are excluded. On the other hand, the red 
lines show the corresponding bounds when only prompt radiation is taken into account. 
The thick black line is obtained when both radiation components are combined for each 
energy regime separately. All fluxes are calculated within our patch S. As obvious from 
this plot ICS radiation bounds dominate at dark matter masses above a few 100 GeV. 




Figure 8. Constraints on decaying dark matter for the decay channel \ — > 
decoded into its different components. The thick solid line shows the overall bounds 
on mass and lifetime, cf. also Fig. 9. Green lines represent the constraint coming from 
the e^-response function for ICS emission alone, whereas red lines are based on the 
prompt photon spectrum alone. Each of the eight lines corresponds to one of the 
observed 7— ray energy ranges as denoted in the caption of Fig. 5. 



In Fig. 9 we show our results for bounds on the four different decay channels into 
H + fi + fi~ , yU + /i~, r + r~ and bb as examples with different amount of prompt and ICS 
radiation. + The first three decay modes can well fit the PAMELA/Fermi positron and 
electron data if the positron excess is interpreted in terms of decaying dark matter, 
and the preferred mass and lifetime regions are indicated by the blue blobs and red 
crosses.* We used the Pythia package [67] to derive the electron, positron and 7-ray 
decay spectra. 

In these plots the dashed-dotted (dotted) line shows the bounds obtained from 
ICS (prompt) radiation in our patch S alone, the thick solid line shows the bounds 
obtained when prompt and ICS radiation are combined. Furthermore, the bounds can 
be strengthened to the yellow region when the foreground model LI is subtracted from 
the data. 

It turns out that for decays into fi + n~ pairs and four-body decay into yU + yU~/i + yU— , 

+ The four-body decay into muons proceeds via two intermediate neutral scalar particles with masses 
of 1 GeV. 

* The shown regions should be understood as typical masses and lifetimes that well fit the data. We 
performed a \ 2 test, keeping the electron background as a freely adjustable powerlaw with spectral 
index between -3.3 and -3.0, whereas the positron background is kept fixed as the one from "Model 
0" in Ref. [32]. We include in the fitting procedure only the PAMELA data on the positron fraction 
above 10 GeV, as they should be less affected by solar modulation, and the Fermi and HESS e ± -data, 
for which we added the corresponding systematic and statistical errors in square. Since systematic 
errors are correlated, the obtained % 2 are relatively small, but in view of the large uncertainties in the 
electron background a more detailed fit is not reasonable. The blue regions correspond to x 2 /dof = 1 
and x 2 /dof = 0.75, outside these regions the fits to HESS and Fermi LAT become problematic. 
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Figure 9. Bounds on different decay channels in the mass vs. lifetime plane. Regions 
below the thick solid line are excluded by combined ICS and prompt radiation in the 
LI propagation model, whereas parameter space below the dashed- dotted (dashed) line 
is excluded due to ICS (prompt) radiation alone. The ICS constraints shown with the 
dashed-dotted lines are calculated from the response functions shown in Fig. 5. The 
constraints can be strengthened to the yellow light shaded region if the predictions 
of Model LI for the galactic diffuse astrophysical foreground is subtracted. The blue 
blobs and red crosses (which are taken from Ref. [21]) show the parameters that well fit 
electron + positron fluxes observed by Fermi LAT and HESS and the positron fraction 
observed by PAMELA as described in the text. 

the strongest constraints typically come from ICS rather than from the prompt radiation 
and the constraints could be improved by more than a factor of 2 for small masses and 
by a few 10% for large masses after removal of the 7-ray emission from conventional 
astrophysical sources. In the case of decay into t + t~ and bb the prompt radiation alone 
already provides strong constraints, which can again be improved by subtracting galactic 
foreground as for the case of decay into muons. 

Note that our patch S is optimized for ICS radiation. Prompt radiation from dark 
matter decay in general dominates at the galactic pole regions, as discussed above (in 
the actual data, this behavior is disturbed at high gamma-ray energies because of the 
large contamination of the data with isotropic cosmic-ray background). Following the 
slicing of the sky as proposed in Ref. [63], we can find for the highest energy bin a patch 
that actually increases our corresponding final state radiation bounds by around 70%. (j 

(j Using this adaptively determined patch, which is located at 10° < b < 20° and < I < 10° and has 
only a few number counts, still does not allow to raise the bounds as high as shown in Ref. [63]. The 
difference might originate in the smaller energy bins used in [63] , and the inclusion of data above 300 
GeV. 
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For comparison, we also show with the dashed lines in Fig. 9 the bounds that can 
be obtained by comparing the sum of extragalactic ICS radiation, extragalactic prompt 
radiation and the maximal isotropic part of the halo prompt radiation (which is identical 
to the flux from the Galactic anti-center) with the preliminary results for the isotropic 
extragalactic gamma-ray flux as presented in Ref. [62]. Comparing these bounds, which 
already rely on a foreground subtraction, in case of decay into muons with the ones 
obtained from patch S after foreground subtraction shows that they are subdominant 
and become only relevant at very high masses. Our bounds are somewhat weaker than 
the ones found in Ref. [66], which is due to our inclusion of absorption effects and our 
more conservative treatment of extragalactic ICS radiation. 

5. Conclusions 

In this study we have calculated the contribution to the 7-ray fluxes from decaying 
dark matter particles, including the inverse Compton photons resulting from energetic 
electrons and positrons through scattering with low energy target photons in addition to 
the bremsstrahlung emissions. We constructed e ± -response functions based on the full- 
sky 7— ray observations by Fermi LAT which can be applied to constrain any decaying 
dark matter model by convolving it with the specific decay spectrum into electrons 
and positrons. We also studied the dependence of the e ± -response functions on both 
the set of propagation parameters and halo profiles and find that the most important 
uncertainty comes from the height of the diffusion zone. 

We applied the response functions to different decaying dark matter models 
with leptonic final states, including also the effects of prompt radiation, which can 
significantly increase the 7— ray emissions in the dark halo. Moreover, we demonstrated 
how the constraints can be further improved by subtracting astrophysical contributions 
to the observed 7— ray flux. 

Under the conservative assumption of a propagation model with the height of the 
diffusion zone around 4 kpc, based on the Fermi LAT data, we can severely constrain 
but not exclude models with dark matter decay into t + t~ that can explain the positron 
excess observed by PAMELA. Moreover, we find that analogous models with two- and 
four-body decay channels into /x ± s remain essentially unconstrained by current obser- 
vations. When our reference foreground model is subtracted the lower bounds on the 
lifetime in general increase by 0(1) factors for dark matter masses below 1 TeV, and 
by 10-60% for masses above 1 TeV, which is however not enough to exclude the above 
channels in the parameter regime relevant for PAMELA. The bounds might improve by 
0(1) factors when data with better background rejection is used. For comparison we 
also calculated conservative bounds from the isotropic extragalactic 7-ray background 
as inferred from the Fermi LAT data, finding again that even the decay into t + t~ can- 
not be excluded in this way. 

In this work we demonstrated the use of e ± -response functions in constraining 
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dark matter models. We plan to update the response functions as data with improved 
rejection of background becomes available. 
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